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^ I The condition of purity of states for a damped harmonic oscihator is con- 

^ • sidered in the framework of Lindblad theory for open quantum systems. For a 

^^ , special choice of the environment coefficients, the correlated coherent states are 

^vi I shown to be the only states which remain pure all the time during the evolution 

of the considered system. These states are also the most stable under evolution 

in the presence of the environment. 
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1 Introduction 

In the last two decades, more and more interest has arisen about the search for a 



^ ! consistent description of open quantum systems [1-5] (for a recent review see ref. [6]). 

O^' Dissipation in an open system results from microscopic reversible interactions between 

the observable system and the environment. Because dissipative processes imply ir- 
r> I reversibility and, therefore, a preferred direction in time, it is generally thought that 

quantum dynamical semigroups are the basic tools to introduce dissipation in quan- 
tum mechanics. In the Markov approximation the most general form of the generators 
of such semigroups was given by Lindblad [7]. This formalism has been studied for 
the case of damped harmonic oscillators [6, 8-12] and applied to various physical phe- 
nomena, for instance, to the damping of collective modes in deep inelastic collisions 
in nuclear physics [13]. A phase space representation for the open quantum systems 
within the Lindblad theory was given in [14, 15]. 

In the present study we are also concerned with the observable system of a har- 
monic oscillator which interacts with an environment. We discuss under what condi- 
tions the open system can be described by a quantum mechanical pure state. In Sec. 
2 we present the generalized uncertainty relations and the correlated coherent states, 
first introduced in [16], which minimize these relations. The Lindblad theory for open 
quantum systems is considered in Sec. 3. Then in Sec. 4, for the one-dimensional 



harmonic oscillator as an open system, we show for a special choice of the diffusion 
coefficients that the correlated coherent states, taken as initial states, remain pure for 
all time during the evolution. In some other simple models of the damped harmonic 
oscillator in the framework of quantum statistical theory [17, 18], it was shown that 
the pure Glauber coherent states remain as those during the evolution and in all other 
cases, the oscillator immediately evolves into mixtures. In this respect we generalize 
this result and also our previous result from [10] as well as the results of other authors 
[19], obtained by using different methods. In Sec. 5 we introduce the linear entropy 
and present its role for the description of the decoherence phenomenon. We claim that 
the correlated coherent states are the most stable under evolution in the presence of 
the environment and make the connection with the work done in this field by other 
authors [20-24]. Finally, concluding remarks are given in Sec. 6. 

2 Generalized uncertainty relations 

In the following we denote by a^^ the dispersion of the operator A, i.e. (T^a =< 
A^ > — < A >^, where < A >= a a = Tr(pyl),Tr/) = 1 and p is the statistical 
operator (density matrix). By aAB = 1/2 < AB + BA > — <A><B>we 
denote the correlation of the operators A and B. Schrodinger [25] and Robertson [26] 
proved for any Hermitian operators A and B and for pure quantum states the following 
generalized uncertainty relation : 

(^AA(rBB -(^\b>^\< [AB] > [^ (1) 

For the particular case of the operators of the coordinate q and momentum p the 
uncertainty relation (1) takes the form 

(Jpp(Tqq - o-lg > — . (2) 

This result was generalized for arbitrary operators (in general non-Hermitian) and for 
the most general case of mixed states in [16]. The inequality (2) can also be represented 
in the following form [16]: 

CTppaqq > ^ (3) 



where 
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is the correlation coefficient. The equahty in the relation (2) is realized for a special 
class of pure states, called correlated coherent states [16], which are represented by 
Gaussian wave packets in the coordinate representation. These minimizing states, 
which generalize the Glauber coherent states, are eigenstates of an operator of the 
form [16]: 

1 r ir 1] 

with real parameters r and rj, \r\ < 1,7] = y/o\^. Their normalized eigenfunctions, the 
correlated coherent states, have the form [16]: 

with a a complex number. If we set r = and r] = {Ti/2mujY^'^, where m and uj are the 
mass and respectively the frequency of the harmonic oscillator, the states (6) become 
the usual Glauber coherent states. In Wigner representation, the states (6) have the 
form [16]: 

1 2rf in — a Y 2r 

This is the classical normal distribution to give dispersion 
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and the correlation coefficient r. The Gaussian distribution (7) is the only positive 
Wigner distribution for a pure state [27]. All other Wigner functions that describe 
pure states necessarily take on negative values for some values of p, q. 

In the case of the relation (1) the equality is generally obtained only for pure 
states [16]. For any density matrix in the coordinate representation (normalized to 
unity) the following relation must be fulfilled: 

7 = Trp2 < 1. (9) 

The quantity 7 characterizes the degree of purity of the state. For pure states 7 = 1, 
for highly mixed states 7 <^ 1 and for weekly mixed states 1 — 7 ^ 1. 

The Wigner function may be expressed as the Fourier transform of the off-diagonal 
matrix elements of the density operator in the coordinate representation: 

W{p, q,t) = \jdy<q- y\p\q + y> e^'^^y'\ (10) 

Tin J 



Then < x\f)\y > can be obtained from the inverse Fourier transform of the Wigner 
function: 

/t X -\- V 

dpexp{-p{x-y))W{p,^—). (11) 

In the language of the Wigner function the condition (9) has the form: 

7 = 27r^ / W^{p, q)dpdq < 1. (12) 

Let us consider the most general mixed squeezed states described by the Wigner 
function of the generic Gaussian form with five real parameters: 

1 1 

W^(P, q) = t; — 7= exp{- — [(7pp(g - cxqf + agg{p - apf - 2apg{q - aq){p - ap)]}, (13) 
Zn^/a Za 

where a is the determinant of the dispersion (correlation) matrix M, 

a = detM = (Tpp(Tgg - a^^ (14) 

and 

M= (""PP ^P^y (15) 

The Gaussian Wigner functions of this form correspond to the so-called quasi-free 
states on the C*-algebra of the canonical commutation relations, which is the most 
natural framework for a unified treatment of quantum and thermal fluctuations [28]. 
For Gaussian states of the form (13) the coefficient of purity 7 is given by 

Therefore, the dispersion matrix has to satisfy the Schrodinger-Robertson uncertainty 
relation (2) 

This inequality must be fulfilled actually for any states, not only Gaussian. Any 
Gaussian pure state minimizes the relation (17). Here, a is also the Wigner function 
area - a measure of the phase space area in which the Gaussian density matrix is 
localized. For a > h /A the function (13) corresponds to mixed quantum states, while 
in the case of the equality a = h /A it takes the form (7) corresponding to pure 
correlated coherent states. 



The degree of the purity of a state can also be characterized by other quanti- 
ties besides 7. The most usual one is the quantum entropy (we put the Boltzmann's 
constant A; = 1): 

5 = -Tr(plnp) = - <lnp> . (18) 

For quantum pure states the entropy is identically equal to zero. It was shown [12, 29] 
that for Gaussian states with the Wigner functions (13) the entropy can be expressed 
through a only: 

5= (z/ + l)ln(z/ + l)-z/lnz/, u=^^--. (19) 

ill Zi 

3 Quantum Markovian master equation for damped 
harmonic oscillator 

The simplest dynamics for an open system which describes an irreversible process 
is a semigroup of transformations introducing a preferred direction in time [2, 3, 7]. 
In Lindblad's axiomatic formalism of introducing dissipation in quantum mechanics, 
the usual von Neumann-Liouville equation ruling the time evolution of closed quan- 
tum systems is replaced by the following Markovian master equation for the density 
operator p(t) in the Schrodinger picture [7]: 

^ = mm. (20) 

Here, $4 denotes the dynamical semigroup describing the irreversible time evolution of 
the open system in the Schrodinger representation and L is the infinitesimal generator 
of the dynamical semigroup $4. Using the Lindblad theorem [7] which gives the most 
general form of a bounded, completely dissipative generator L, we obtain the explicit 
form of the most general quantum mechanical master equation of Markovian type: 

^ = -\\iirpit)\ + ^E([^.pW>^/] + \%rp{t)v}\)- (21) 

Here H is the Hamiltonian operator of the system and V^-, K are bounded operators 
on the Hilbert space 7i of the Hamiltonian and they model the environment. We make 
the basic assumption that the general form (21) of the master equation with a bounded 
generator is also valid for an unbounded generator. In the case of an exactly solvable 
model for the damped harmonic oscillator, we define the possible two operators Vi and 
V2, which are linear in p and g, as follows [6, 8, 9]: 

V, = a,v^h,q, J = 1,2, (22) 



with aj, bj complex numbers. The harmonic oscillator Hamiltonian H is chosen of the 
general form 

H = Ho + f^iqp + pq), Ho = ^p' + ^e. (23) 

With these choices and with the notations 

h h h 

Dqq = i^Y. \(^3?^Dpp = 7: Y. l^jf^^pq = Dqp = -^tRc ^ a*bj,X = -Im ^ a*bj,{2A) 

i=l,2 ^ j=l,2 ^ j=l,2 j=l,2 

where a* and b* denote the complex conjugate of Qj, and bj, respectively, the master 
equation (21) takes the following form [6, 9]: 

do t ^ i i 

^ = —^[Ho,p] - ^{^ + f^)[q, pp + pp] + ^(A-/i)[p,pg + gp] 

-^[g, [q,p]] - ^[p, [p,p]] + |p([g, [p,p]] + [p, [q,p]]). (25) 

Here the quantum mechanical diffusion coefficients Dpp, Dqq, Dpq and the friction con- 
stant A satisfy the following fundamental constraints [6, 9]: Dpp > 0,Dqq > and 

DppDqq -Dl^>^. (26) 

The relation (26) is a necessary condition for the generalized uncertainty inequality 
(2) to be fulfilled. 

The semigroup method is valid for the weak-coupling regime, with the damping 
A typically obeying the inequality A ^ cjq, where cjq is the lowest frequency typical of 
reversible motions. 

The necessary and sufficient condition for L to be translationally invariant is 
/i = A [6, 8, 9]. Translation invariance means that [p, L{p)\ = L{[p, p]). In the following 
general values for A and p will be considered. 

By using the fact that the linear positive mapping defined by A ^ Tr(pA) is 
completely positive, in [6, 9] the following inequality was obtained: 

DppCTqqit) + Dqqappit) - 2Dpq(Tpg(t) > —— . (27) 

We have found in [10] that this inequality, which must be valid for all values of t is 
equivalent with the generalized uncertainty inequality (2) at any time t, 

^qqit)app{t) - alq{t) > -, (28) 



if the initial values (7qg{0),app{0) and o"pg(0) for t = satisfy this inequality. If the 
initial state is the ground state of the harmonic oscillator, then 

^,M = 1^^, ^Pp(O) = ^-, ap,(0) = 0. (29) 

By using the complete positivity property of the dynamical semigroup <l>i, it was 
shown in [9] that the relation 

IV($,(p)5:v/v;) = ^Tr($,(p)V^t)Tr($,(^)p..) (gg) 

3 3 

represents the necessary and sufficient condition for p(t) = $t(/3) to be a pure state for 
all times t > 0. This equality is a generalization of the pure state condition [30-32] to 
all Markovian master equations (21). If p^(t) = p(t) for all t > 0, there exists a wave 
function ip &7i which satisfies the nonlinear Schrodinger-type equation 

th^ = H'm, (31) 

with the non- Hermit ian Hamiltonian 



H' = H+tj2< m,v,H{t) >v,~-< m,T.y3V3m > -^ E^/^r (32) 
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For environment operators Vj of the form (22), the pure state condition (30) takes the 
following form [9], corresponding to equality in the relation (27): 

Dpp(Jqq{t) + Dqqapp{t) - 2Dpq(Tpq{t) = —— (33) 

and the Hamiltonian (32) becomes 

H' = H + X{ap{t)q - aq{t)p) + ^[Xfi^ - Dqq{{p - ^^(t))^ + app{t)) - Dpp{{q - aq{t)f 

+(rqq{t)) + Dpqiip - ap{t)){q - aq{t)) + (g - a,(t))(p - ap{t)) + 2apq{t))]. (34) 

From the master equation (25) we obtain the following equations of motion for 
the expectation values and variances of the coordinate and momentum: 

^ = -(A-p)a,(t) + la,(t), (35) 

at m 

dap{t) 



dt 



moj (Jq{i) — (A + n)ap{t) (36) 



and 



^^ = _2(A - f,)a,,it) + ^a,,{t) + 2D,„ (37) 



-^ = -2(A + ^)app(t) - 2mcc;Vpg(t) + 2Dpp, (38) 

'^^^ = -muj\g{t) + -app{t) - 2Xap,{t) + 2D„. (39) 

LLL lit 

In the under damped case {oo > fi) considered in this paper, with the notation Q^ = 
iJ^ — /i^, we obtain [6, 9]: 

ag(t) = e~^*((cosl]t + ^sinl]t)(Tg(0) + — - sinf]tc7p(0)), (40) 

i L I I tu U 
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ap{t) = e"^*( — sinfita,(0) + (cosl^t- ^ sinfit)ap(O)) (41) 

and o"g(oo) = (7p(oo) = 0. It is convenient to consider the vectors 

X{t)=\ a„,{t)/mu I (42) 

and 

(2mujDqq 
2Dpp/muj I . (43) 

2Dpq 

With these notations the solutions for the variances can be written in the form [6, 9]: 
X{t) = (Te^*T)(X(0) - X(cx))) + X(oo), (44) 

where the matrices T and K are given by 

, / fi + iQ fi — iQ 2uj \ 

^-iVt fi + iVt 2uj , (45) 

\ —uo —uo — 2/i/ 



2iVt 



-2{\-iVl) 

K=\ -2{\ + m) I (46) 

-2A 
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and 

X(oo) = -{TK-^T)D. (47) 

The formula (47) is remarkable because it gives a very simple connection between 
the asymptotic values (t -^ oo) of aqq{t),a.pp{t),a.pq(t) and the diffusion coefficients 

^qqi ^'ppi PI ■ 

-,,(oo) = 2(^^)2a(a'+c.^-^^) (("^")'('^(^ + ^) + "')^- 

+cu2Dpp + 2mcu2(A + /i)Dpg), (48) 



(y'nrAOO] 



^ " \2. .2n , /o\/\ ..\ , .2\ n o™. .2/ 



^^'^^' 2A(A2 + CU2 - ;i2) 



{{muoyuj'Dqq + (2A(A - yu) + cj^)Z^pp - 2mcj^(A - ;u)Dpg),(49) 



f^P9(oo) = 2m\i\^ + uj^- 2) ^"^^ + /i)(m^)2D^, + (A - ^)Dpp + 2m(A2 - ^i^)Dpq).{hQ) 

These relations show that the asymptotic values (Tgg(oo), (Tpp(oo), (jpq{oo) do not depend 
on the initial values o-qg(O), o-pp(O), o'pg(O). In the considered underdamped case we have 

-2At / ^11 ^12 ^13 

^21 &22 &23 I , (51) 

V ^31 ^32 ^33 , 



Te^^T 



2^2 

where 6jj, i,j = 1,2,3 are time- dependent oscillating functions given by (3.78) in [9]. 

4 Purity of states 

We will be interested to find the Gaussian states which remain pure during the evo- 
lution of the system for all times t. We start by considering the pure state condition 
(33) and the generalized uncertainty relation (2) which transforms into the following 
minimum uncertainty equality for pure states: 

By eliminating (jpp between the equalities (33) and (52), like in [33], we obtain: 



/ /,x _ Dpq(^pq(t) + 4^ ^ .2 , ^PP^QQ ~ ^pq r/ /,x 4^ ^Dpq 

y^qqV^J n ) ' n2 [{^pqV') 



^PP ^pp ^pp^qq ^pq 

1 2 -Ppp-Pgg--Pp9- 4^ ^ n n 1-n i^-X) 

+ 4 m n _ n2 ^2 ^pp^qq\-^- K.^^} 
^ K-^pp-^qq ^pq) 



Since the opening coefficients satisfy the inequahty (26), from the relation (53) we 
obtain the following relations which have to be fulfilled at any moment of time: 

D,pD,,-Dl^ = ^, (54) 



Dpp(Tgg{t) - Dp^dpy^t) — = 0, (55) 



ap,{t){DppD,, - Dl^) - ^Dp, = 0. (56) 

From the relations (52) and (54) - (56) it follows that the pure states remain pure for 
all times only if their variances are constant in time and have the form: 

M^) = ^, M^) = ^, cr,,(t) = :^. (57) 

If these relations are fulfilled, then the inequalities (27) and (28) are both equivalent 
to (26), including also the corresponding equalities (33), (52) and (54). From Eq. 
(44) it follows that the variances remain constant and do not depend on time only 
if X{0) = X(oo), which means o-pp{0) = (Tpp(oo), <yqq{0) = 0-^^(00), o"pg(0) = 0-^^(00). 
Using the asymptotic values (48) - (50) of the variances and the relations (57), we 
obtain the following expressions of the diffusion coefficients which assure that the 
initial pure states remain pure for any t : 

Formulas (58) are generalized Einstein relations and represent examples of quantum 
fluctuation-dissipation relations, connecting the diffusion with both Planck's constant 
and damping constant. With the coefficients (58), the variances (48) - (50) become: 



2mQ' """ ^?r' '^^^""21^' 



'^11 o™n' '^PP on ' ^P9 on- y^^) 



Therefore, the quantity a (see Eq. (14)) is equal to its minimum possible value h /4, 
according to the generalized uncertainty relation (17). Then the corresponding state 
described by a Gaussian Wigner function is a pure quantum state, namely a correlated 
coherent state [16] with the correlation coefficient (4) r = —^ijuj. Given Opp^ Oq^ and 
(Tpq, there exists one and only one such a state minimizing the uncertainty o [34]. A 
particular case of our result (corresponding to A = /i and Dpq = 0) was obtained by 
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Halliwell and Zoupas by using the quantum state diffusion method [19]. We remark 
that the minimization of the quantity (14) is equivalent, by virtue of the relation 
(19), to the minimization of the entropy S. As we have mentioned before, we have 
considered general coefficients fi and A and in this respect our expressions for the 
diffusion coefficients and variances generalize also the ones obtained by Dekker and 
Valsakumar [33] and Dodonov and Man'ko [35], who used models where /i = A was 
chosen. If yU = 0, we get from (58) Dpg = 0. This case, which was considered in [10], 
where we have obtained a density operator describing a pure state for any t, is also a 
particular case of the present results which give the most general Gaussian pure state 
which remains pure for any t. For /i = 0, the expressions (59) become 

which are the values of variances (29) for the ground state of the harmonic oscillator 
and the correlation coefficient (4) takes the value r = 0, corresponding to the case of 
usual coherent states. 

The Lindblad equation or its equivalent Fokker-Planck equation for the Wigner 
function with the diffusion coefficients (58) can be used only in the underdamped case, 
when uj > fi. Indeed, for the coefficients (58) the fundamental constraint (26) implies 
that m^(cj^ — ii^)D'l > /i^A^/4, which is satisfied only ii cu > fi. It can be shown that 
there exist diffusion coefficients which satisfy the condition (54) and make sense for 
CO < fi, but in this overdamped case we have always a > Ti /A and the state of the 
oscillator cannot be pure for any diffusion coefficients [35]. 

The fluctuation energy of the open harmonic oscillator is 

1 1 

E{t) = —crpp{t) + -muj^agg{t) + fiapg{t). (61) 

If the state remains pure in time, then the variances are given by (57) and the fluctu- 
ation energy is also constant in time and is given by 

11 1 

E = ji^Dpp + -mU^Dgg + fiDpg). (62) 

Minimizing this expression with the condition (54), we obtain just the diffusion coeffi- 
cients (58) and Emm = fiQ/2. Therefore, the conservation of the purity of state implies 
that the fluctuation energy of the system has all the time the minimum possible value 
Emin- Moreover, we can consider that the case when the diffusion coefficients satisfy 
the equality (54) corresponds to a zero temperature of the environment (bath). Then 
the influence on the oscillator is minimal. 

11 



If we choose the coefficients of the form (58), then the equation for the density 
operator can be represented in the form (21) with only one operator V, which up to a 
phase factor can be written in the form: 



with [V, 1/t] = 2h\. 

The correlated coherent states (6) with nonvanishing momentum average, can also 
be written in the form: 

^(a^) = (t; )^ exp[-- (1 - -^ap,)ix - a,y + ^apx] (64) 

and the most general form of Gaussian density matrices compatible with the general- 
ized uncertainty relation (2) is the following: 



IN / 1 \i r 1 ,^ + y 

< x\p\y >= (ttt^)' ^^pHtt^I^:^ ^q 



i2 



' 27ro-<j, 2aqq 



^2 



-T^^i^pp - -^)i^ - yf + F^(^^ ^9)(^ -y) + T^p{x - y)]. (65) 

2ft (Jqq haqq 2 Tl 

For these matrices we can verify that Trp^ = Tijl^ and they correspond to the 
correlated coherent states (64) if a^p, Oqq and u^q in (65) satisfy the equality in (2). 

In [11, 14, 15] we have transformed the master equation for the density operator 
into the following Fokker-Planck equation satisfied by the Wigner distribution function 

dW V dW n dW d d 

+B„_ + C^^_ + i,„_. (66) 

For an initial Gaussian Wigner function the solution of Eq. (66) is 

X exp{- — [app(g - aq{t)y + aqq{p - (Jp{t)f - 2apq{q - (rq{t)){p - (Tp{t))]}. (67) 

We see that the initial Wigner function remains Gaussian and therefore the property of 
positivity is preserved in time. Consider the harmonic oscillator initially in a correlated 
coherent state of the form (64), with the corresponding Wigner function (7). For our 
environment described by the diffusion coefficients (58), the solution for the Wigner 
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function at time t is given by (67), where (Tq{t) and o'p(t) are given respectively by (40) 
and (41) and the variances by (59). Using either (11) and (67) or (65), we get for the 
density matrix the following time evolution: 



,exp-^- -T^-^,(^^^' 



< x\m\y >= (^)^ exph — (:^ - a,it)f 



mil, ,r, imu,x + y , ,, , . i , ., ^^ ,^r 

'x - yf - -^{^ - ^,(t))(x - y) + -ap{t){x - y)]. (6^ 



Ah' ''^ h ' 2 '' ''' "' n 

In the long time limit (Tg(t) = 0, o"p(t) = and we have 

< x|p(cx))|i/ >= (^)^ exp{-^[f](x2 + 1/2) + ,^(a;2 _ y2^^^^ (gg) 

TiTi 2h 

The corresponding Wigner function has the form 

Wooip, 9) = ^ """"P'^M] ^m ^ "^^'^' ^ ^'"^^^^' ^'^^^ 

We see that the time evolution of the initial correlated coherent state of the 
damped harmonic oscillator is given by a Gaussian density matrix with variances 
constant in time. According to known general results [22, 23], the initial Gaussian 
density matrix remains Gaussian centered around the classical path. So, the correlated 
coherent state remains a correlated coherent state and crq{t) and cXpit) give the average 
time dependent location of the system along its trajectory in phase space. 

5 Entropy and decoherence 

Besides the von Neumann entropy S (18), there is another quantity which can measure 
the degree of mixing or purity of quantum states. It is the linear entropy Si defined as 

Si = Tt{p-p^) = 1-Ttp\ (71) 

For pure states Si = and for a statistical mixture Si > 0. 

As it is well-known, the increasing of the linear entropy Si (as well as of the 
von Neumann entropy S) due to the interaction with the environment is associated 
with the decoherence phenomenon (loss of quantum coherence), given by the diffusion 
process [20, 21]. Dissipation increases the entropy and the pure states are converted 
into mixed states. The rate of entropy production is given by 

Slip) = -2Tr(p/)) = -2Tr(pL(p)), (72) 

where L is the evolution operator. According to Zurek's theory, the maximally predic- 
tive states are the pure states which minimize the entropy production in time. These 
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states remain least affected by the openness of tlie system and form a "preferred set 
of states" in the Hilbert space of the system, known as the "pointer basis". Their 
evolution is predictible with the principle of least possible entropy production. 
Using (25), in our model the rate of entropy production (72) is given by: 

4 
Slit) = -^[DppTTip^^ - MM) 

+^<?qTr(pV - PPPP) - DpgTiip'^iqp + pq) - 2Mpp) - — Tr(/32)]. (73) 

When the state remains approximately pure (p^ ~ p), we obtain: 

Slit) = j^iDppCTggit) + D,,a„it) - 2Dp,ap,it) - — ) > 0, (74) 

according to (27). We see that Siit) = when the condition (33) of purity for any time 
t is fulfilled. The entropy production Si is also equal to at t = if the initial state 
is a pure state. But we have shown before that the only initial states which remain 
pure for any t are the correlated coherent states. Therefore, we can state that in 
the Lindblad theory for the open quantum harmonic oscillator the correlated coherent 
states, which are generalized coherent states, are the maximally predictive states. Our 
result generalizes the previous results which assert that for many models of quantum 
Brownian motion in the high temperature limit the usual coherent states correspond to 
minimal entropy production and, therefore, they are the maximally predictive states. 
In our model the coherent states can be obtained as a particular case of the correlated 
coherent states by taking /i = 0, so that the correlation coefficient r = (see Eq. (4)). 
Paz, Habib and Zurek [20, 21] considered the harmonic oscillator undergoing quan- 
tum Brownian motion in the Caldeira-Leggett model and concluded that the minimiz- 
ing states which are the initial states generating the least amount of von Neumann or 
linear entropy and, therefore, the most predictible or stable ones under evolution in the 
presence of an environment are the ordinary coherent states. Using an information- 
theoretic measure of uncertainty for quantum systems, Anderson and Halliwell showed 
in [22] that the minimizing states are more general Gaussian states. Anastopoulos and 
Halliwell [23] offered an alternative characterization of these states by noting that these 
states minimize the generalized uncertainty relation. According to this assertion, we 
can say that in our model the correlated coherent states are the most stable states 
which minimize the generalized uncertainty relation (2). Our result confirms the one 
of [23], but the model used in [23] is different, namely the open quantum system con- 
sists of a particle moving in a harmonic oscillator potential and is linearly coupled to 
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an environment consisting of a bath of harmonic oscillators in a thermal state. At 
the same time we remind that the Caldeira-Leggett model considered in [20, 21] vio- 
lates the positivity of the density operator at short time scales [36, 37], whereas in the 
Lindblad model considered here the property of positivity is automatically fulfilled. 

The rate of predictibility loss, measured by the rate of linear entropy increase, for 
a damped harmonic oscillator is also calculated in the framework of Lindblad theory 
in Ref. [38]. The initial states which minimize the predictability loss are identified as 
quasi-free states with a symmetry dictated by the environment diffusion coefficients. 
For an isotropic diffusion in phase space, the coherent states or mixtures of coherent 
states are selected as the most stable ones. 

In order to generalize the results of Zurek and collaborators, the entropy pro- 
duction was considered by Gallis [24] within the Lindblad theory of open quantum 
systems, treating environment effects perturbatively. Gallis considered the particular 
case with Dpq = and found out that the squeezed states emerge as the most stable 
states for intermediate times compared to the dynamical time scales. The amount of 
squeezing decreases with time, so that the coherent states are most stable for large 
time scales. For Dpg 7^ we have generalized the result of Gallis and established that 
the correlated coherent states are the most stable under the evolution in the presence 
of an environment. 

6 Concluding remarks 

Recently there is a revival of interest in quantum Brownian motion as a paradigm of 
quantum open systems. The possibility of preparing systems in macroscopic quantum 
states led to the problems of dissipation in tunneling and of loss of quantum coherence 
(decoherence) . These problems are intimately related to the issue of quantum to 
classical transition and all of them point the necessity of a better understanding of 
open quantum systems. The Lindblad theory provides a selfconsistent treatment of 
damping as a general extension of quantum mechanics to open systems and gives 
the possibility to extend the model of quantum Brownian motion. In the present 
paper we have studied the one-dimensional harmonic oscillator with dissipation within 
the framework of this theory. We have shown that the only states which stay pure 
during the evolution in time of the system are the correlated coherent states under 
the condition of a special choice of the environment coefficients. These states are also 
connected with the decoherence phenomenon and they are the most stable under the 
evolution in the presence of the environment. In a future work in the framework of the 
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Lindblad theory we plan to discuss the connection between uncertainty, decoherence 
and correlations of open quantum systems with their environment in more details. 
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